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ABSTRACT 

We report on the highly variable Si IV and C IV broad absorption lines in SDSS 
J113831.4+351725.2 across four observational epochs. Using the Si IV doublet compo¬ 
nents, we find that the blue component is usually saturated and non-black, with the ra¬ 
tio of optical depths between the two components rarely being 2:1. This indicates that 
these absorbers do not fully cover the line-of-sight and thus a simple apparent optical 
depth model is insufficient when measuring the true opacity of the absorbers. Tests 
with inhomogeneous (power-law) and pure-partial coverage (step-function) models of 
the absorbing Si IV optical depth predict the most un-blended doublet’s component 
profiles equally well. However, when testing with Gaussian-fitted doublet components 
to all Si IV absorbers and averaging the total absorption predicted in each doublet, the 
upper limit of the power law index is mostly unconstrained. This leads us to favour 
pure partial coverage as a more accurate measure of the true optical depth than the 
inhomogeneous power law model. 

The pure-partial coverage model indicates no significant change in covering frac¬ 
tion across the epochs, with changes in the incident ionizing flux on the absorbing 
gas instead being favoured as the variability mechanism. This is supported by (a) the 
coordinated behaviour of the absorption troughs, (b) the behaviour of the continuum 
at the blue end of the spectrum and (c) the consistency of photoionization simula¬ 
tions of ionic column density dependencies on ionization parameter with the observed 
variations. Evidence from the simulations together with the C IV absorption profile in¬ 
dicates that the absorber lies outside the broad line region, though the precise distance 
and kinetic luminosity are not well constrained. 
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1 INTRODUCTION 

Blue-shifted quasar absorption lines, which indicate the 
presence of large amounts of outflowing material from the 
central engine, have become an increasingly important topic 
over the past two decades. This followed the realisation of 
the potential influence these flows ca n exert on the host 
g alaxy by r estricting black hole growth (|Springel et al.|[2005i : 
iKind l20ld l and contributing to the evolution of galac¬ 
tic bulges llMaeorrian et alJ Il998l : Isilk fe Rees] 1 1 9981 ) . This 
feedback process is also supported by evidence that out¬ 
flows could extend hun dreds of pars ecs into the surround¬ 
ing interstellar medium (lBarlowlll994h , de positing significan t 
amounts of kinetic energy in the process (lAr av et a l. l2013f) . 

Absorption features observed in the ultraviolet (UV) 
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portion of the quasar’s rest-frame spectrum can be blue- 
shifted by as much as 0.1c, where c is the speed of light 
in a vacuum, from the corresponding emission line centre. 
These absorption lines are divided into categories according 
to their velocity width, the widest being broad absorption 
lines (BALs) extending over at least 2000 km s^ 1 where flux 
is les s than 90 per cent of the continuum ( Wevman n et al.l 
EaaS)- Quasars hosting at least one BAL are known as broad 
absorption line quasars (BALQSOs). Absorption can also be 
seen as narrow absorption lines (NALs), which extend over 
a few hundred km s -1 and mini-BALs, which are of in ter¬ 
medi ate width between NALs and BALs dHamann fe Sabral 
l2004h . 

The population of BALQSOs, of which SDSS 
J113831.4+351725.2 (hereafter SDSS J1138+3517) is a 
member, account for approximately 10 per cent to 20 per 
cent of the total quasar population llReichard et al.l 120031 : 
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iKnigge et al.ll2008l : IScaringi et al.ll2009h . However SDSS tar¬ 
get select ion effects may m ake this value as high as 41 
per cent (lAllen et al.l l201l] '). If BAL observability is gov¬ 
erned by line of sight orientation, then large-scale out¬ 
flows may occur in the v ast majority of quasar s, even 
where no BALs are visible (jSchmidt feHineall999h . Com¬ 
mon BAL subtypes include High-Ionization BALQSOs (Hi- 
BALs) (~85 per cent) whose spectra only show broad ab¬ 
sorption due to high ionization transitions such as C IV 
AA1548,1550 Si iv AA1394,1403 and N v A A1239,1243 
llSpravberrv fe Foltz! 1 19921 ; iReichard et al.l l2003h . Low Ion¬ 
ization BALQSOs (LoBALs), in addition to high-ionization 
BALs, also exhibit broad absorption resulting from low ion¬ 
ization species such as Mg II, A1 ill and on rare occasions 
Fe ii (known as FeLoBALs). 

Ionizing continuum changes a re known to drive 
broad emission line (BEL) variability llPeterson et al .1 l~i~998l : 

lYanden Berk et aI]|2004[Wilhite et alJhood h however there 


is as yet no consensus behind the dominant mechanism 
involved in absorption line variablity. Several large sta¬ 
tistical studies of BAL variability have been published 


in recent years (Lundgren et al. 

2007 

: Gibson et al. 

2008 

20101; Capellupo et al. 201ll. 2012 

2013 

; Filiz Ak et al. 

2013 

Wildv et al. 20141). These studies refer to coordinated vari- 


ability across widely separated absorption troughs as evi¬ 
dence favouring changes in inp ut ionizing flux as th e main 
mechanism behind variability dFiliz Ak et al.l l2013l l , while 
variability over narrow velocity ranges within troughs is at- 
tributed to absorbing gas moving across the line of sight 


( Gibson et al. 2008; Capellu 

po et al. 2012|). 

As in Wildv et al. (2011 

), we use a novel method to re- 


construct the un-absorbed spectral profile (I Allen et al.ll20lih 
against which the absorption can be measured. Absorption 
is identified based on this emission profile as it contains the 
total input spectral flux (continuum+emission line) enter¬ 
ing the absorbing gas at each wavelength bin. We use this 
method to estimate the column density at each epoch for 
the Si IV and C iv ions. This has advantages over purely 
continuum-based absorber strength measurements as it al¬ 
lows absorption depth to be accurately identified in spectral 
regions where the absorption is overlapping with an emis¬ 
sion line profile, allowing absorption features to be studied 
to arbitrarily low (less negative) velocities. 

1.1 The Unusual BAL behaviour of SDSS 
J1138+3517 

The quasar examined in this paper, SDSS J1138+3517, has 
a redshift of z=2.122 (jHewett fe W ildll 201 oil and was investi¬ 
gated as part of the multi-BALQSO variab ilit y stu dy under¬ 
taken in the rest-frame UV bv lWildv et al.l (l2014T l. Using a 
sample of 50 quasars, Si iv A1400 and C iv A1549 BALs were 
identified and statistics gathered on their equivalent width 
(EW) changes over 2 epochs separated by quasar rest-frame 
timescales ranging from approximately 1 to 4 years. Of the 
50 BALQSOs, SDSS J1138+3517 showed the largest abso¬ 
lute change in EW in both Si IV and C iv absorption as well 
as a large fractional EW change (| A EW / {EW) |), occurring 
over a rest-frame time interval of ~1 year. A third epoch 
indicated another dramatic variation in the BALs, while a 
fourth epoch showed relatively little change from the third 
except for significant variability in one C IV absorber, with 


these subsequent observations also being made in successive 
1-year rest-frame time intervals. 

In iFiliz Ak et al.1 (120131 1 . the variability of 428 C iv 
BALs in 291 quasars was found to produce a A EW de¬ 
pendence on rest-frame time interval which is well predicted 
by a random-walk model (see Fig. 29 of that paper). Us¬ 
ing the same BAL identification scheme as that study, the 
variation acr oss the two observational epochs examined in 
IWildv et al.l (1201414 of the most variable C iv BAL region 
in SDSS J1138+3517 (spanning —23 400 km s” 1 to —7500 
km s -1 ) is AEW=— 13.8+2.9 A. This value lies well outside 
the random-walk model prediction at the epoch separation 
time-interval of 360 days in the quasar rest-frame, suggesting 
that this object may exhibit a variability mechanism which 
is rare in the BALQSO population. 

In lFiliz Ak et all (120141 4. which uses the same BAL iden¬ 
tification scheme as above, the variability of C IV BALs 
was investigated over timescales ranging from ~1 to a few 
years, including those in quasars additionally containing 
Si iv BALs of overlapping velocity. Again examining the 
most variable C iv BAL in SDSS J1138+3517 identified us¬ 
ing their definition, the value of A EW, along with a value of 
AEW/ (EW)=— 1.3+0.3 would place it among the top few 
most variable C IV BALs in their sample of 454 (includes 
C IV BALs both unaccompanied by Si IV absorption and 
overlapping in velocity with a Si iv BAL), highlighting the 
extreme variab il ity of this so urce. In addition, an investi¬ 
gation bv lGibson et al.1 (|201Q|4 included observations of two 
quasars over a similar rest-frame time interval of between 0.5 
and 2 years, neither of which showed variability as dramatic 
as in SDSS J1138+3517. 

Answering the question of what drives the strong BAL 
variability observed in this quasar is the principal objective 
of this investigation. Observations of the source are detailed 
in Section 2. In Section 3, we quantify Si iv and C iv ionic 
column densities using different assumptions regarding the 
line-of-sight geometry, including models where there is in¬ 
complete coverage of the emission region by a homogeneous 
absorber. This is aided by the use of Gaussian profiles to fit 
the Si iv absorption components, which allows calculations 
to be performed based on the changes in these components 
even where they are not fully un-blended. In Section 4 we 
compute a grid of photoionization models, allowing exam¬ 
ination of the theoretical variation of ionic column density 
over a range of ionization parameter at given hydrogen num¬ 
ber densities and column densities. These grid models are 
compared to the best results for the ionic column densities 
from Section 3, allowing us to test the hypothesis of changes 
in ionizing flux being the driver of BAL variability. In Sec¬ 
tion 5 we discuss the results from Sections 3 and 4 and their 
implications for the variability mechanism. We also attempt 
to constrain the mass outflow rate and hence kinetic lumi¬ 
nosity. 

This paper assumes a flat ACDM cosmology with 
Ho=67 km s -1 Mpc -1 , Gm=0.315 and Ba=0.685, as 
reported from the latest results of the Planck satellite 
(IPlanck Collaboration et alj|20l3 ). 
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2 OBSERVATIONS 

2.1 The Sloan Digital Sky Survey 

The epoch 1 observation was obtained from DR6 of the 
Sloan Digital Sky Survey (SDSS). Beginning in 2000, the 
SDSS had i maged approximately 10 000 deg 2 of the sky as of 
June 2010 (ISchneider et al.lfeoid 'l. using the 2.5m Telescope 
at th e Apache Point Observatory, New Mexico, USA dYorkl 
l200dh Imaging is carried out using a CCD cam era that op¬ 
erates in drift-scanning mode rtGunn et al.lll99§ ). Five broad 
band filters, u,g,r,i and z, normalised to the AB system, are 
used, covering a wavelength range of 3900 to 9100 A with an 
approximate spectral r esolution of A/AA=2000 at 5000 A 
llStoughton et all 120021 '). Further information on the SDSS 
data reduction process is given in iLupton et alJ ll200ll ') and 
IStoughton et all (I2002T ). 


2.2 The William Herschel Telescope 

The William Herschel Telescope (WHT) ISIS (Intermedi¬ 
ate dispersion Spectrograph and Imaging System) double¬ 
armed (red and blue) spectrograph was used for longslit 
spectroscopy of the target during semester A 2008, 2011 
and 2014 corresponding to epochs 2,3 and 4 respectively, 
with epochs 3 and 4 being observed as part of the ING ser- 
vice programme an d epoch 2 taken from the sample used by 
ICottis et al.l ()2010i 'l. The observations provide spectral cov¬ 
erage from 3000 to 10 000 A. For epoch 2 the 5700 dichroic 
was used, providing spectra from the two arms with over¬ 
lap in the range 5400 to 5700 A. The R158B and R158R 
gratings, giving a nominal spectral resolution of 1.6A per 
pixel in the blue and 1.8A per pixel in the red, were also 
utilised for epoch 2. Epoch 3 and 4 observations employed 
the 5300 dichroic, giving an overlap range of 5200-5500 A, 
while the R300B and R316R gratings were used for the blue 
and red observations, providing spectral resolutions of 0.86 
A per pixel and 0.93 A per pixel respectively. In the dichroic 
overlap regions, an error weighted mean was used to com¬ 
bine the spectra from the two arms for each observation. 
For all WHT observations, a slit width of 1.5” was cho¬ 
sen to match typical seeing conditions and give reasonable 
throughput without compromising the spectral resolution. 
Multiple exposures were taken (in order to remove cosmic 
ray hits) bracketed by arc-lamp exposures for wavelength 
calibration. Standard star spectra were taken on the observ¬ 
ing nights to allow flux calibration of the targets. Correction 
of CCD images using bias and flat field exposures was per¬ 
formed within the iraf ccdproc routine. Spectra were ex¬ 
tracted, along with sky background removal, using the apall 
task within the iraf longslit package. This used an extrac¬ 
tion slit width of 10 pixels and optimal weighting, with 2 
pixels per wavelength bin. Wavelength and flux calibration 
were applied to the extracted spectra within this same pack¬ 
age. The spectra were placed on an absolute flux scale using 
spectra of a photometric standard observed at similar air- 
mass. No attempt was made, such as through the use of grey 
shifts, to adjust the flux scale on the calibrated spectra to¬ 
wards levels which more closely matched across the epochs. 
Table [T] summarises the SDSS J1138+3517 observations. 
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3 ANALYSIS 

3.1 Spectral Properties 

The observed spectra at each epoch correspond to rest-frame 
UV wavelengths which include the region containing Si IV 
and C IV emission and the corresponding blueshifted absorp¬ 
tion lines resulting from those same atomic transitions in 
outflowing gas. Since the spectral resolution differs between 
spectra, in order to accurately identify changes in absorp¬ 
tion between each epoch the SDSS and WHT epoch 3 and 
epoch 4 spectra were convolved with Gaussians of appro¬ 
priate full-width half maxima to approximate the resolution 
of the WHT epoch 2 observation (using skyline widths as 
a guideline), followed by resampling of all 4 spectra onto a 
wavelength grid of 3.7 A. This closely matches the bin width 
of the epoch 2 observation. Small differences in wavelength 
calibration were removed by aligning narrow emission and 
absorption features in the WHT spectra with the same lines 
in the SDSS spectrum. The SDSS and both WHT spectra 
were corrected for the effects of_Milky Way Galactic Ex¬ 
tinction using the method of Cardelli et al.l ( 1989l l with Ay 
values taken from ISchlegel et alJ (jl998h . 

Using the reconstructi o n o f the unabsorbed SDSS 
spectrum fr om Wildv et al.l (120141 ') originally developed in 
lAllen et al.l d201lf> . appropriate reconstructions were also 
generated for the three WHT observations. This was 
achieved as follows: (i) A power-law continuum was fitted 
to each spectrum us ing line-free continuum bands as out¬ 
lined in Wildv et al.l (l2014li . (ii) The SDSS continuum was 
subtracted from the SDSS reconstruction, (iii) The WHT 
spectra had their respective continua subtracted, (iv) In the 
continuum subtracted SDSS reconstruction, the Si iv A1400, 
C iv A1549, O 1 A1304 and C 11 A1334 emission line bound¬ 
aries were identified, (v) The emission lines from part (iv) 
were scaled to match the red side (unabsorbed) of the emis¬ 
sion lines in each of the WHT observations, and (vi) The 
continuum was added to the scaled emission lines for each 
of the WHT observations, creating appropriate reconstruc¬ 
tions for epochs 2,3 and 4. 

For the SDSS reconstruction used in the subsequent 
analysis, the procedure was carried out to scale the Si iv 
emission line, as the original reconstruction underestimates 
the Si iv emission (t he reconstructions were optimised for 
C iv as described in [Alien et al.l (1201 ill 1. The region used 
for calculating the appropriate C iv line scaling is restricted 
to the interval between line centre and an observed wave¬ 
length of 4866 A. This is to avoid the He Il+Om absorption 
and emission complex immediately longward of C iv. The 
Si IV emission line was shifted by two wavelength bins blue- 
ward in the epoch 2 reconstruction relative to the SDSS 
reconstruction to match the shifted peak of the emission 
line in the epoch 2 spectrum. Based on the synthetic BAL 
based error estimation method described in Section 3 of 
[Wildv et all (l2014f) . the fractional error on the reconstruc¬ 
tion for a quasar of this r-band S/N and redshift is estimated 
to be 5 per cent of the flux value at each resolution element 
for all calculations described in this paper. The observed 
spectra for all 4 epochs (solid black lines), as well as their fi¬ 
nal reconstructions (solid red lines), are illustrated in Fig.Q] 

The dramatic variability in both ions’ absorbers across 
the first three epochs is clear when the spectra are nor¬ 
malised to the reconstructions. The transitions of interest 
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Table 1. Details of observations obtained of SDSS J1138+3517 


Epoch 

Observation Date At qres t 
(days) 

Type 

Exp. time 

00 

Grating 

Pixel Resolution (5000 A) 
(km s -1 ) 

Mean Airmass 

1 

03 May 2005 

0 

SDSS 

6x2400 

Red & Blue 

150 

1.01 

2 

27 May 2008 

359 

WHT 

3x1200 

R158B & R158R 

108 

1.14 

3 

29 May 2011 

710 

WHT 

2x1800 

R300B & R316R 

56 

1.80 

4 

09 May 2014 

1053 

WHT 

2x1800 

R300B & R316R 

56 

1.06 


A tqrest is quasar rest-frame time interval since Epoch 1 


Rest-frame wavelength (A) 

1300 1400 1500 1600 



Observed wavelength (A) 


Figure 1. Spectra for epochs 1 to 4 of SDSS J1138+3517, with 
top panel to bottom panel in order of observation date starting 
with the earliest. The observed spectrum is in black while the re¬ 
construction of the unabsorbed spectrum is in red. Vertical dotted 
lines indicate the laboratory wavelengths of the Si IV and C IV 
emission lines. The vertical dashed line indicates the maximum 
wavelength used in calculating the C IV emission scaling, emission 
longward of this point is not accurately reconstructed. 


are both doublets, with rest-frame laboratory wavelengths 
of 1393.76 and 1402.77 A contributing to Si iv A1400 and 
correspondingly 1548.20 and 1550.77 A for C iv A1549. Rela¬ 
tive to the laboratory-frame rest-wavelength of the red com¬ 
ponent of each doublet, absorption is significant between 0 
and 13 000 km s -1 for Si iv and between 0 and ^—20 000 
km s" 1 for C iv, where a negative value indicates outflow¬ 
ing material (blueshifted). As can be seen in Fig. [2] the 
C IV troughs are deeper than the Si iv troughs at similar 
velocities. The deepest troughs show the least variability 
and are located at the lowest (least negative) velocities, as 



Velocity (km s ! ) 


Figure 2. Si IV (upper panel) and C IV (lower panel) absorption 
regions in red component velocity space. Epoch 1 (SDSS) spectra 
are in black, epoch 2 (1st WHT) spectra are in red, epoch 3 (2nd 
WHT) spectra are in green and epoch 4 (3rd WHT) spectra are in 
blue. Spectra are normalised to the un-absorbed reconstruction. 
The Si IV Gaussian components (described in Section 3.3) are 
numerically labelled in order of increasing outflow velocity 


has been reported in previou s BAL studies (jLundgren et al.1 
120071 : ICapellupo et al.1 [201 lh . 

Due to its relative lack of variation from the previous 
observation, especially in the Si iv absorption lines where 
identifiable doublets are present, epoch 4 is left out of the 
analysis in subsequent parts of Section 3 and Section 4. 
Epoch 4 is instead discussed further in Section 5.3. 
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3.2 Lower Limits for Outflow Column Densities 
using Direct Integration 

Minimum values for both the Si IV and C IV column den¬ 
sities can be estimated by assuming unsaturated absorbers 
completely cover the emitting line+continuum region along 
our line of sigh t with a constant optical depth across the tan¬ 
gential plane dSavage fc Sembachlll99lT ). The C IV and Si IV 
transitions have a doublet structure with known oscillator 
strengths. The total column density of the ion contributing 
to a given doublet component may be calculated according 
to 

N ion = [ r{y)dv (1) 

7re 2 /A J 

where Ni on is the ionic column density, v is the velocity 
relative to the laboratory rest-frame wavelength, m e is the 
electron mass, c is the speed of light, e is the elementary 
charge, / is the oscillator strength, A is the laboratory rest- 
frame wavelength and t(v) is the optical depth at a given 
velocity. At a redshift of z = 2.122, the doublet separation at 
zero velocity is approximately 28.1 A for Si IV and 8.0 A for 
C IV. As the C IV separation is only slightly greater than two 
wavelength bins, it is practically impossible, given their large 
intrinsic widths, to identify individual C IV components. 
Though some Si IV components are well-separated, many 
of the doublets are blended to various extents with neigh¬ 
bouring absorption features. However, the optical depth of 
overlapping components will be the sum of the individual 
components and the optical depth of the red component is 
known to be half that of the blue component for both ions, 
given the ratio of their oscillator strengths. Therefore, by 
integrating over all of the absorption for a given transition 
using the value of / corresponding to the red component, a 
minimum value of Ni on can be found by taking 1/3 of the 
calculated value. The velocity ranges integrated over and 
resulting estimated ionic column densities are shown in Ta¬ 
ble [2] 

Since the C IV doublet is unresolved, this method is 
the only means of estimating a lower limit for the column 
density of this ion. For Si IV, methods involving modelling 
the components need to be checked for consistency with the 
lower limit derived here. 


3.3 Gaussian Components of the Si IV Outflow 

The Si IV doublet velocity separation is ~1920 km s -1 , 
allowing some individual components to be identified. As¬ 
suming narrow components of quasar absor ption lines fol¬ 
low an approximately Gaussian profile, e.g. lHamann et al.l 
Hoil] ). model profiles can be constructed for the doublet 
lines. Spectral fitting is performed using the specfit package 
within iraf, requiring three free parameters for each Gaus¬ 
sian component, namely the wavelength at line centre, line 
full width half maximum (FWHM) and line EW. Appropri¬ 
ate restrictions are applied to these parameters, i.e. the EW 
ratio between the blue and red components of the doublet 
must be between 2:1 and 1:1, the difference in velocities at 
line centre must not differ from the expected doublet sepa¬ 
ration by more than one velocity bin, and the widths must 
be the same. By using an initial ’guess’ for these parameters’ 
values, spectral fitting is performed using chi-square minimi- 



-8000 -6000 -4000 -2000 0 


Velocity (km s _1 ) 

Figure 3. Epoch 1 (top panel), epoch 2 (middle panel) and 
epoch 3 (lower panel) models of Si IV doublet absorption lines 
in red component velocity space (related doublet components are 
shown in the same colour). The black dotted line describes the 
total model absorption profile, while the red dotted line near the 
bottom of each panel indicates the difference between the model 
profile and the observed profile. 


sation. In total, five Si IV doublets are identified across the 
3 epochs, spanning the range — lOOOO^v^O km s -1 . These 
are shown in Fig. [3] 

Comparing the Gaussian fit to the almost un-blended 
red component of the lowest velocity doublet in epoch 2 
(shown in red in Fig. [BJ to the data points at these veloc¬ 
ities gives a fractional error in the fit of 0.033. To be con¬ 
servative a larger error of 5 per cent on each point of all 
Gaussian profiles is assumed for calculations involving these 
model profiles. While the EW and occasionally the FWHM 
change from one epoch to the next, there is no significant 
change in the line centre velocity and hence no evidence for 
acceleration of the outflowing material. This lack of accelera- 
tion has been noted in pr evious studies (iHamann et al.|[2008l : 
iRodrfguez Hidalgo et al.l l2011). The best-fit parameters ob¬ 
tained by specfit for each doublet component are provided 
in Table [3] 

Upon examination of an individual component, the fol¬ 
lowing equation can be used to estimate the column density 
giving rise to the doublet to which it belongs: 


N = 


m e cb 

y/ne 2 f\ 


to 


( 2 ) 


where tq is the peak optical depth and b is the Doppler pa- 
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Table 2. C IV and Si IV velocity limits and ionic column densities for each epoch using direct integration of the absorption profile. 


Transition 

Velocity Limits 

Epoch 1 Ionic Column Density 

Epoch 2 Ionic Column Density 

Epoch 3 Ionic Column Density 


(km s _1 ) 

(xl0 14 cm -2 ) 

(xl0 14 cm -2 ) 

(xl0 14 cm —2 ) 

Si iv A1400 

-13100 

24.7i2.66 

7.21i2.37 

11.7i2.41 

C iv A1549 

-19 700 

190T15.2 

95.8i8.64 

151ill.3 


t Integration spans zero velocity to velocity limit 


Table 3. List of components with parameters for Gaussian model profiles. Doublets are listed in order of increasingly negative velocity 
(1 is lowest, 5 is highest). Doublet velocity is measured from the red component rest wavelength. 


Absorber 

Velocity-^ 
(km s -1 ) 

Epoch 1 

Epoch 2 

Epoch 3 



FWHMt 

EW+ 

FWHMt 

EWt 

FWHMt 

EWt 



(km s _1 ) 

(A) 

(km s -1 ) 

(A) 

(km s _1 ) 

(A) 

1 (Red) 

-2300 

900 

3.18i0.09 

900 

2.03i0.07 

800 

2.69i0.04 

1 (Blue) 



3.18i0.16 


2.80i0.07 


3.21i0.05 

2 (Red) 

-3300 

800 

1.65A0.12 

600 

0.32i0.07 

900 

1.OliO.06 

2 (Blue) 



1.76A0.08 


0.60i0.08 


1.62i0.07 

3 (Red) 

-4300 

1100 

2.62i0.24 

500 

0.19i0.08 

900 

0.60i0.11 

3 (Blue) 



2.92i0.18 


0.19i0.07 


0.72i0.07 

4 (Red) 

-4800 

1000 

2 .21i0.22 

500 

0.16i0.08 

700 

0.48i0.07 

4 (Blue) 



2.88i0.18 


0.21i0.08 


0.97i0.06 

5 (Red) 

-5900 

700 

0.77i0.15 

700 

0.15i0.08 

700 

0.09i0.06 

5 (Blue) 



1.02i0.13 


0.15A0.11 


0.16i0.06 


TEW in quasar rest-frame 
iError ~200 km s _1 


rameter (related to line width by 6 = %/2 x FWHM/2.355). 
This assumes a Voigt profile for the absorbers, which tends 
towards the Gaussian case in the limit of low optical depth. 
Red components are less likely to be affected by satura¬ 
tion since they are weaker than their corresponding blue 
components. As a result, column density measurements 
are taken from red doublet lines. As in Section 3.2, es¬ 
timates of the absorbing column density are lower limits 
as the absorber may not cover 100 per cent of the emis¬ 
sion region. Summing the values obtained for the five dou¬ 
blets listed in Table 0 gives estimates of the total col¬ 
umn densities as Nsi iv,i=33.5il.49, Nsi rv,2=7.79A0.50 
and Nsi iv,3=14.7i0.47 in units of 10 14 cm -2 for the first, 
second and third epochs respectively. As expected, these are 
larger than the limits in Table [2] since blue:red doublet com¬ 
ponent ratios are less than 2:1. 


3.4 A Pure Partial Coverage Model 

Many absorption line systems have been found to be best de¬ 
scribed by a situation in which material of a certain optical 
depth covers a fraction of t he emission source, leaving the 
rest o f the source uncovered (lArav et al.ll 19991 : Ide Kool et al.l 
1200 j) otherwise known as a pure partial coverage (PPC) 
model. Under these conditions the apparent optical depth 
r = ln(l/I res ), where 7 res is the residual intensity (nor¬ 
malised to the unabsorbed continuum+line emission flux), 


will underestimate the true optical depth of the absorbing 
gas and consequently will underestimate the column density. 
The observed flux I app is related to the covering fraction C 
by Ia PP = Clout + (1 — C)Io, where Iq is the unabsorbed 
emission flux and lout is the true output flux from the ab¬ 
sorber. Thus the covering fraction C as a function of velocity 
is given by 


, . = [I (V) - 2 I r (v) + 1 
1 ' I b (v) - 2 Ir ( V ) + 1 


(3) 


where l r ( v ) and I(, ( v ) are the apparent residual fluxes of the 
red and blue components respectively at a given velocity v. 
Since the true optical depth r = ln(Io/I 0 ut), the following 
is also true 


r 


■ 


/ Ir — Ib 
l 1 -Ir 


) 


(4) 


where r is the optical depth of the red component at a spe¬ 
cific velocity. Since (7^1, the values of l r and I5 can be shown 
to be constrained to the range I r ^ It, ^ 7 2 . Assuming this 
condition is met, it is then straightforward to calculate the 
column density using Equation [T) 

Evidence for pure partial covering can be obtained by 
comparing the residual flux velocity dependence of a (al¬ 
most) saturated absorption line to the profile of 1 — C. Simi¬ 
lar profiles would suggest the shape of the absorption profile 
is determined by the velocity-dependent covering fraction of 
an opaque absorber rather than intrinsic differences in opti- 
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Velocity (km s ') 


Figure 4. The absorber 1 profiles of the red and blue compo¬ 
nents interpolated onto the red component’s velocity grid. The 
solid red and blue lines represent the Gaussian models, the dot¬ 
ted lines represent the observed fluxes, the green short-dashed 
line indicates the 1 — C values at these points and the magenta 
long-dashed line represents the e~ T values of the red component. 
Appropriately coloured triangles and squares represent the ob¬ 
served points in velocity space. While the 1 — C profile is a good 
match for the blue absorber, the e~ T profile does not match the 
red profile. 


cal depth. An additional test is to compare the profile of e~ T 
to the red component of an absorption doublet. If the pro¬ 
files do not match, the absorption line profile cannot be due 
to an absorber which is homogeneous and completely cov¬ 
ers the emission region. In order to carry out these tests, a 
doublet with unblended components is desirable. The closest 
example available in the dataset is absorber 1 in the epoch 
2 spectrum. From the values of the EWs of the red and blue 
components (see Table [BJ, the blueired ratio is less than 2:1, 
indicating that the blue component is probably saturated. 
After placing the blue component onto the red component 
velocity grid by interpolation, six velocity bins are found 
which follow the general shape of the corresponding Gaus¬ 
sian models for both the red and blue components. Five of 
these points satisfy IrMb^Ir and are illustrated in Fig. [7] 
From the profile traced by 1 — C (Fig. [4] dashed green 
line), it is clear that the saturated blue component could 
have its profile determined by the uncovered fraction. It is 
also obvious that the residual intensity of the red compo¬ 
nent is a poor match to e~ T . Together these provide strong 


evidence that partial coverage is a significant effect in this 
absorber and possibly others in this quasar. 

Given that this is the only doublet where meaningful 
results can be obtained on an individual velocity bin ba¬ 
sis, and even in this case it can only be done close to the 
peak of the absorption components, it is necessary to use 
the Gaussian models to achieve estimates of column density 
across a greater proportion of the Si IV absorption for the 3 
epochs. Individual velocity bins in the Gaussian models have 
no physical realisation, so unlike for the previous example it 
is impossible to construct 1 — C profiles using Equation [3] 
However if for each doublet component, model fluxes are 
averaged over all velocities spanned by the absorption, in¬ 
dividual components can be treated as individual velocity 
bins. This method can subsequently be used to estimate av¬ 
erage values of covering fraction and true optical depth for 
each doublet. Using this true optical depth in Equation [l] 
gives the doublet’s column density along the line of sight 
when Ni 0 „ is multiplied by the covering fraction. These val¬ 
ues are listed in Table [4] Due to the weakness of absorbers 3 
to 5 in epochs 2 and 3 this can only be performed across all 
absorbers in epoch 1, since at small absorber depths opti¬ 
cal depth and covering fraction measurements are unreliable 
due to their sensitive dependence on doublet component ra¬ 
tio. 

3.5 An Inhomogeneous Absorber Model 

An alternative to the PPC model is a situation whereby the 
optical depth varies with spatial location across the line of 
sight. This was originally develo ped using a power-law de¬ 
pendent optical depth model in |dejvool_et^d . d2002l l and 
was subsequently investigated by lArav et alj (l2005i ). who 
found power-law and Gaussian shaped optical depth depen¬ 
dencies on a 1-dimensional spatial location to be inadequate 
in describing O VI absorbers in Mrk 279. Here, we attempt 
a similar investigation using a power-law model of the form 
r(x,A)=T maa ;(A)x a , where r(x,A) is the optical depth at po¬ 
sition x on the surface of the absorber and at wavelength A, 
t max (A) is the maximum optical depth of the absorber and 
a is the power -law i ndex. _ 

Following lArav et alj (120051 1. we make two simplifying 
assumptions to aid calculations without any loss of general¬ 
ity. First, we assume that the surface brightness is uniform 
across the face of the emission source such that S(x,y,A)=l, 
where x and y are spatial positions in the plane of the line 
of sight and S is arbitrarily normalised. We also simplify the 
optical depth spatial dependence by using a 1-dimensional 
model, so that r(x,y,A)=r(x,A). The x values are fixed to 
span an interval O^x^l for simplicity. The observed resid¬ 
ual flux at a given velocity bin will then be found using 
J-res — fo* T ^dx. To assess a broad range of power-law in¬ 
dices, values from 0 to 10 with step size of 0.1 are tested. In 
order to integrate over the x range, it is effectively divided 
into 1000 bins, giving d:r=0.001. 

Similar to the procedure in Section 3.4, we first exam¬ 
ine the plausibility of a power-law model describing the real 
behaviour of the absorption features by attempting to pre¬ 
dict the profile of the blue component of absorber 1 in epoch 
2 after interpolation of the blue absorption line profile into 
the velocity space of the red component. Given a value of a 
to be tested, the process is carried out as follows: (i) Find 
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Table 4. Table of covering fractions and column densities calculated for absorbers 1 to 5. Absorbers 3 to 5 are not listed for epochs 2 
and 3 due to their low strength. 


Absorber 

Epoch 1 

Epoch 2 

Epoch 3 

c 

Nion C 

Nion C 

Nion 


X 10 14 cm —2 

xl0 14 cm -2 

xl0 14 cm -2 


1 

0.225T0.013 

90 0+15.6 
°—14.0 

0.234T0.042 

7 1 7+ 3 -56 
' " L ' —2.87 

0.236T0.020 

19 o+4.40 
±z *°—3.881 

2 

0.187T0.017 

11 o+ 8 - 68 

J " L,U —7.40 

0-318“ 

0 77‘i+ 13 -7 

u. / / o_ 0 773 

0.275T0.165 

2*898±S : °£° 

3 

0.339T0.017 

14 1+4.20 





4 

0.366T0.032 

0 /ic+2.34 
8.4P_ 2 09 





5 

0.122T0.043 

2 Q4 - *" 3 - 61 






Table 5. Power indices calculated by matching the predicted 
blue component of an inhomogeneous power-law to the Gaussian 
model equivalent (’unc’ indicates ’unconstrained’). Several of the 
values have reached the upper limit of 10, necessarily leaving the 
upper bound unconstrained in these cases. 


Absorber 

Power-law index 


Epoch 1 

Epoch 2 

Epoch 3 

1 

10.0±“ n 9 c 

q q +unc 

io.o±“7 

2 

10.0+“7 

1-7 -1 -\-unc 

a r +unc 
D,D -3.6 

3 

10.0±“" c 

c A+unc 

0 -2.9 

1 rj+imc 

1 ' u -0.4 

4 

6A -1.7 

3 . 9^7 

4.6tr 2 c 

5 

/-V rj-\-UnC 

u - ' -0.4 

0 n~\~unc 
S-d-5.1 

io.oj:“” c 


the value of Tmax which gives the value of the residual red 
component flux I r =J Q e~ TmaxX dx at each velocity bin, (ii) 

Predict each /(, value using lb = J 1 e~ 2TmaxX dx, and (iii) 
Compare predicted h values to the observed values and cal¬ 
culate the corresponding y 2 value over the entire velocity 
range. This procedure is repeated for all values of a, the 
value adopted being the one which minimises the y 2 calcula¬ 
tion in step (iii). This provides a best value of a=3.3^Q'j, the 
corresponding predicted blue profile is shown in Fig. |3] Find¬ 
ing Ni on requires knowledge of the average opti cal depth t, 
where f = f r(x)dx. As in Fig. 2 and Fig. 3 oflA rav et al.l 
(l2005h . the dependence of the doublet component residual 
intensities on f as well as comparisons with complete homo¬ 
geneous coverage and PPC models are shown in Fig. [6] 

From Fig. [5] it is clear that, like the PPC model, an 
inhomogeneous (power-law) absorber model is a good pre¬ 
dictor of the behaviour of the blue component in this par¬ 
ticular doublet. As in Section 3.4, this is extended to the 
Gaussian models with doublet components treated as sin¬ 
gle points with flux averaged over all velocity bins showing 
absorption. Interestingly, this leaves the upper bound on a 
unconstrained for all absorption doublets except absorber 4 
at epoch 1 (see Table [5]). 

Given that the best-fit power-law index of the Gaus¬ 
sian model for absorber 1 at epoch 2 is not consistent with 
the measured value of a from the six velocity bins in the 
observed spectrum, this analysis is considered inconclusive. 
The constraints on the strongest doublet (absorber 1) sug¬ 
gest that this feature has a comparatively large value of a. 
At large values of a the dependence of optical depth on posi- 



Velocity (km s ') 


Figure 5. Inhomogeneous absorber model (red dotted line) of the 
blue doublet component for a power-law index a=3.3. The solid 
blue line is the blue component Gaussian model and the dotted 
blue line is the observed spectrum in velocity range of the blue 
component, both of which have been interpolated onto the red 
component velocity grid. Power-law model error bars are based 
on the errors of a. Blue triangles and red squares indicate the 
points on the observed spectrum and inhomogeneous absorber 
model respectively. 


tion, x , approximates a step function, and thus cannot easily 
be distinguished from a PPC model. 
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Figure 6. Behaviour of doublet components with respect to aver¬ 
age optical depth for the power-law model adopted (a=3.3). The 
residual intensity of the blue component I), (blue line) is found 
by doubling the optical depth of the simulated red component I r 
(red line). The if value (magenta line) indicates what the blue 
component would have been for a homogeneous absorber com¬ 
pletely covering the emission source. The green line is (1-C) for a 
PPC model given the values of I r and It, at a specific f. The ra¬ 
tio °f indicated by the black line, where rppc is the optical 

depth of the PPC model predicting the (1 — C) curve, shows the 
divergence of the power-law model from the PPC model. 

4 PHOTOIONIZATION SIMULATIONS 
4.1 Cloudy Setup 

Modelling the column densities and ionization fractions of 
ions in quasar outflows by specifying the shape of the ioniz¬ 
ing continuum and various gas (column) densities can pro¬ 
vide insight into the physical location of the outflowing 
plasma and the changing nature of the relationship between 
the ionizing photon flux and ionization state, especially if 
ionization changes are the principal driver of absorption 
line variability. Photoionization models in this investigation 
are performed using CLOUDY (iFerland et al.lfl998l ). version 
cl3.02. 

The CLOUDY code combines these input parameters 
with a model of the input spectral energy distribution (SED) 
to calculate the total ionizing photon flux. It models the in¬ 
put SED using input parameters for temperature ( T ) and 
power-law indices determining continuum behaviour of the 
form F v oc v a , including optical to near-UV spectral index 
( a uv ) for optical wavelengths above 2500 A, UV to soft x-ray 


spectral index ( a ox ) between 2500 A and 2 keV, and x-ray 
spectral index (a x ) from 2 keV to 100 keV. The a uv value 
used is 0, which is w ithin the range found by other studies, 
for example iNatali et alj (Il998l ) found an average value of 
a uv =— 0.33±0.59. We ad opt a v a lue of a ox = —1.98 calcu¬ 
lated using the method of lWilkes et al.l (j 19941 ) in Equation [5] 

a M =-1.53-0.11 xlog(j^l) (5) 

where L 0 is the specific intensity at 2500 A in the quasar 
rest-frame in units of erg s -1 Hz -1 . This is towards t he more 
negative end of the range found in lGrupe et al.l (|2010| ). how¬ 
ever tests with values of a ox between a more typical value 
of —1.6 and the adopted value show no significant difference 
in Si IV and C IV column densities, so any discrepancies 
between our adopted value and the true va lue shou ld not 
drama tically alter our results. According to Zdziarski et al. 
d 19961) the value of a x is typically between —0.8 and —1.0, so 
a value of —1.0 is adopted here for simplicity. It is possible 
to estimate T from the ratio of the bolometric to Eddington 
lu minosity and the mas s of the black hole using the method 
of lBonning et al.l ll2007l) : 

T ma , = 10 5 ' 56 A/-^i^ (6) 

LEdd 

where T max is the maximum accretion disk temperature, 
Mg is the black hole mass in units of 10 8 Mq, Lt, 0 i is the 
bolometric luminosity and Lndd is the Eddington luminosity. 
The value of LEdd can be found from the black hole mass 
using Lsdd = 10 46 1 Mg. Using estimates of Mg and Lbol from 
iShen et ahl (l201ll ) gives a value of T max r^ 240 000 K, with 
the peak flux of the big blue bump occurring at 2.1xl0 16 
Hz (~140 A or 87 eV). 

A grid of models containing column density data for 
Si IV and C IV is generated using the input continuum, 
with grid-points at hydrogen column densities log (Nh / 
cm~ 2 )=21, 22 and 23, hydrogen number densities log (uh / 
cnr 3 )=5, 7 and 9 and log U values spanning —5.0 to 3.0 in 
intervals of 0.2 dex. These hydrogen number densities span 
typical narrow line region to broad line region (BLR) densi¬ 
ties, the maximum density being close to the critical density 
of C ill]. The density is unlikely to be higher than this as this 
would suggest the absorber originates from within the BLR, 
conflicting with the fact that individual absorption compo¬ 
nents are relatively narrow in comparison to the broad emis¬ 
sion lines. This is also apparent from the need for the deepest 
troughs to be absorbing both the BEL and the continuum. 
The ionization parameter range encompasses all reasonable 
values given the existence of the absorbers, while the hy¬ 
drogen column density range was chosen to span a range 
above the minimum ionic column density for C IV, Nc iv, 
and below a density at which Thomson scattering becomes 
significant. This rules out column densities of log (Nh / 
cm~ 2 )^20 as these generate Nc iv values below the mini¬ 
mum value calculated from direct integration (Fig. Q. 

4.2 Application of Cloudy Models to Estimated 
Column Densities 

The CLOUDY models can be applied to results from the pre¬ 
vious sections in this paper to provide parameters governing 
the properties of the quasar continuum source which affect 
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Figure 7. Predicted C IV column densities as a function of hy¬ 
drogen number density for given hydrogen column densities. Solid 
squares indicate the CLOUDY input hydrogen densities, which are 
linked by solid lines of constant hydrogen column density. The 
horizontal dashed line gives the minimum C IV column density 
predicted by the SDSS observation, which shows the strongest 
absorption of all epochs. It is clear that log (Nh / cm -2 )^20 is 
ruled out by this limit. 


the outflowing gas. Initially it is necessary to find the correct 
normalisation of the continuum SED described in Section 
4.1, which can be used to find the total ionizing luminos¬ 
ity Li on and ionizing photon emission rate Q(H ) and hence 
the distance from the continuum source to the inner edge of 
the outflowing gas (facing the source) and the mass outflow 
rate. The correct ionising SED is found by scaling the out¬ 
put to the rest-frame corrected flux of the observed SDSS 
spectrum. Output parameters from the simulation include 
the total ionizing photon flux (f>(H ) which can be used to 
find the true values of the ionising luminosity, the rate 
of emission of Lyman continuum photons, Q(H ) and the 
ionising source distance to the absorbing gas, R. 

The calculated value of R will vary depending on 
the grid parameters (U,Nh and ng); however since the 
same SED is used in each case, values of ionizing lu¬ 
minosity and ionizing photon production rate are found 
which are applied universally; they are Li on =4.17x 10 46 
erg s _1 and Q(H)=9.80x 10 56 s -1 . The elemental abun¬ 
dances us ed by cloudy are a s olar composition de- 
rived from iGrevesse fc Sauvall (1 1 9981 ) , iHolweaen d200lTl , and 
lAllende Prieto et al. (2001, 2002th 

In order to apply the grid of models calculated in Sec- 


Table 6. Total column densities and their limits for Si IV at each 
epoch. 


Epoch 

Asi iv 

(xlO 14 cm -2 ) 

JVsi iv lower limit 
(xlO 14 cm -2 ) 

iVgi iv upper limit 
(xlO 14 cm -2 ) 

1 

59.8 

33.5 

94.2 

2 

9.18 

7.79 

27.5 

3 

18.3 

14.7 

27.9 


tion 4.1, limits for the column densities of Si IV and C IV 
must be found. Since there are no resolvable doublets avail¬ 
able for C IV, only the lower limits found from direct inte¬ 
gration (see Section 3.2) can be used for this ion. For Si IV, 
the strength of the doublets at epoch 1 means total column 
densities can be estimated using the PPC model for all 5 
absorbers as listed in Table [4] For epochs 2 and 3, PPC 
can be used for absorbers 1 and 2 which are the predomi¬ 
nant contributors to the total line of sight column density 
at these epochs, plus a contribution from absorbers 3 to 5 
using the peak optical depth (POD) method outlined in Sec¬ 
tion 3.3 modified to take account of partial coverage. Here 
we assume that the covering fractions of absorbers 3 to 5 
at epochs 2 and 3 are the same as those at epoch 1. The 
true optical depth can be found once the covering fraction 
is known. Multiplying the original POD-derived Nsi iv by 
the covering fraction and the ratio of the integrated real op¬ 
tical depth over the integrated apparent optical depth span¬ 
ning the velocity range of the red component of the doublet, 
gives an estimate of the true column density. The Nsi iv and 
their upper and lower limits are shown in Table [S] The total 
column density lower limits are assumed to be the totals de¬ 
rived from the POD method (listed in Section 3.3). Upper 
limits on the value of Nsi iv at epoch 2 are comparatively 
large due to the large uncertainty in the covering fraction 
of absorber 2, however the best value lies towards the lower 
end of this interval. 

Using the CLOUDY output at each combination of 
log Nh and log tih, plots can be made of predicted ionic 
column densities across the entire range of ionization pa¬ 
rameter U. The span of log U covers the peak of both Ni on 
for all combinations, resulting in cases where there are two 
possible regions in which the ionic column density values 
can be located. Given that only lower limits of the C iv 
column density can be found, the column densities for this 
ion are estimated by locating the positions in log U space of 
allowed Nsi iv and using the equivalent values for C IV. If 
it is assumed that ionizing continuum changes are the main 
factor driving variability, it is possible to identify the range 
of log U which satisfies the ionic column density limits at 
each epoch. Although there are nine possible combinations 
of log Nh and log tih, only 2 are shown (Figures [8] and [9]). 
It is apparent from these Figures that the lower region of al¬ 
lowed log U values (where the gradient of log Ni on vs. log U 
is positive) is invalid as the allowed range of Nsi iv for epoch 
3 does not permit any Nc iv values above the lower limit. 
This was found to be true for all 9 density/column density 
combinations, giving one continuous span of allowed log U 
values in each case. 
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log U 


Figure 8. Predicted Si IV (solid red lines) and C IV (solid blue 
lines) column densities as a function of ionization parameter at 
log (Nji / cm _2 )=22 and log (n# / cm _3 )=7. Black, red and 
green dashed lines with arrows indicate the lower limits of the 
C IV column density for epochs 1, 2 and 3 respectively. Using the 
same colour scheme, the solid triangles represent the best values 
of the Si IV column densities, while the dashed-dotted lines show 
the confidence intervals on these values with the area enclosed 
being shaded by vertical dotted lines. The thick solid lines at the 
top show the corresponding log U limits. 

4.3 Outflow properties derived from models 

The estimation of mass outflow rates ( M ou t) is based on the 
method used in lCrenshaw et ahl (|2009l ) which relates M ou t 
to other parameters using the following equation 

Mout = 8ttRNh n m pC q \v\ (Moj/r -1 ) (7) 

where /x is the mean atomic mass per proton (/x~ 1 ), m p is 
the proton mass, C q is the covering fraction of the outflow 
as seen from the quasar and v is the mean velocity of the 
outflow. The value of C q is assumed to be ~0.1 given the 
fraction of quasars showing broad absorption lines. It is not 
obvious what the mean outflow velocity of the predominan- 
tely hydrogen outflow is, so we have adopted the centroid 
velocity of the C IV outflow, v~—7600 km s -1 . 

In Table [7] information on mass outflow rate, distance 
to the inner face of the outflowing gas from the ionising 
continuum source and C IV column density is provided at 
the values of U allowed by the limits on Nsi iv and lower 
limits on Nc iv at each grid point. This gives an overview of 
the allowed range in physical state of the outflow, as well as 
possible effects a changing ionization parameter could exert. 



log U 


Figure 9. Predicted Si IV (solid red lines) and C IV (solid blue 
lines) column densities as a function of ionization parameter at 
log (Nji / cm _2 )=23 and log (nj. j / cm _3 )=9. Black, red and 
green dashed lines with arrows indicate the lower limits of the 
C IV column density for epochs 1, 2 and 3 respectively. Using the 
same colour scheme, the solid triangles represent the best values 
of the Si IV column densities, while the dashed-dotted lines show 
the confidence intervals on these values with the area enclosed 
being shaded by vertical dotted lines. The thick solid lines at the 
top show the corresponding log U limits. 


Linear interpolation between the 0.2 dex interval-separated 
log U values is used to estimate values at intervals of 0.02 
dex. 

As expected, the range of log U allowed is strongly de¬ 
pendent on Nh rather than uh, leaving the value of the 
hydrogen number density uncertain, with the caveat that it 
probably does not exceed the BLR density. If ionizing con¬ 
tinuum flux variations are driving the variability, then given 
the values of log U and their possible ranges indicated by 
the third and fourth columns of Table [7] it is possible to es¬ 
timate how the ionization parameter is changing across the 
time scales sampled. Between epochs 1 and 2, log U will be 
changing by increments of between 0.06 dex and 0.76 dex 
with best estimates of ~0.5 to ~0.6 dex. For the epoch 2 to 
epoch 3 interval, log U changes by increments of between 0 
and 0.38 dex, with 0.2 dex being the best estimate. This 
would require a similar change in the ionizing flux reaching 
the outflowing gas. Changes of similar scale to these best 
estimates for the change in log U have been observed in 
the extreme ultraviolet (EUV) of at least one other quasar 
of similar redshift over a comparable rest-frame timescale 
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Table 7. Table of parameters derived from CLOUDY simulations for various hydrogen densities and column densities. 


log N h 

log n H 

log U 

log U range 

log N c IV 

log N c jy range 

log R 

Outflow Rate 

log(cm —2 ) 

log (cm -3 ) 



log(cm —2 ) 

log(cm —2 ) 

log (cm) 

M©/year 

Epoch 1 

21 

5 

-1.20 

-1.36<log 1.06 

17.09 

16.91<log N c rv<17.12 

20.8T0.1 

36.lt®' 4 

21 

7 

-1.24 

-1.38<log t/^-1.06 

17.10 

16.88^1og N c i V ^17.15 

19.8T0.1 

3.6itr 7 t 

21 

9 

-1.18 

-1.34<log C7<-1.02 

17.07 

16.85^1og Nq iv <17.15 

18.8T0.1 

0.36lt°°?4 

t22 

5 

-0.20 

-0.30<log t/^-0.08 

17.90 

17.32^1og N c IV ^18.03 

20.3T0.1 

114^23* 

22 

7 

-0.24 

-0.34^1og tr^-0.10 

17.87 

17.33$log N c iv<18.03 

19.3T0.1 

n-4±l;g 

22 

9 

-0.20 

-0.30<log t/<-0.08 

17.76 

17.22^1og N c i V ^17.94 

18.3T0.1 

1 14+ 0 - 30 

0.23 

t23 

5 

0.80 

0.68<log [7^0.92 

18.79 

17.66sSlog N c i V ^18.98 

19.8T0.1 

36lt®4* 

23 

7 

0.76 

0.66<log U ^0.90 

18.72 

17.73$log N c iv<18.94 

18.8T0.1 

36.lt®' 4 

23 

9 

0.80 

0.70^1og U^O.92 

18.48 

17.48i£log N c iv<18.76 

17.8T0.1 

3-6it8;?l 

Epoch 2 

21 

5 

-0.64 

-0.98^1og U^-0.62 

16.37 

16.31<log N c i V <16.81 

20.5j:°i 

18.1 + 3° 7 6 

21 

7 

-0.66 

-1.00<log U^-0.62 

16.34 

16.29<log Nc iv^16.80 

19.5+°; 2 

1 01+1.05 
-0.37 

21 

9 

-0.62 

-0.94<log U^-0.58 

16.33 

16.28<log N c i V ^16.74 

I8.5±g;f 

n Ifii +0- 106 
U.181_o 037 

+22 

5 

0.26 

-0.02<log RjgO.30 

16.46 

16.40<log N c i V ^17.03 

20.0tJJ;? 

57 4+ 33 - 3 

22 

7 

0.24 

-0.04^1og 0.26 

16.44 

16.42sglog N c iv<17.04 

19 1+ 0 - 2 

7 20~*~ 4 ' 20 
' ,ZU -1.48 

22 

9 

0.28 

-0.02^1og t/^0.30 

16.42 

16.39sSlog N c iv<16.95 

18.0tg' 2 

0 574+0 333 

+23 

5 

1.32 

l.OO^log U 0.36 

16.33 

16.27sSlog N c iv<16.90 

19.5«; 2 

181^?®* 

23 

7 

1.30 

0.98<log U 0.34 

16.32 

16.27i£log N c i V ^17.00 

l8.5+°; 3 

18.lt3 8 7 ° 

23 

9 

1.34 

1.00<log U 0.38 

16.29 

16.23^1og N c iv<16.81 

17.5 t° 0 i 

1 01 +1-06 

1-8 -0.37 

Epoch 3 

21 

5 

-0.86 

-0.98^1og C/<-0.82 

16.63 

16.57^1og N c i V <16.81 

20.6T0.1 

22.8^® 

21 

7 

-0.88 

-1.00<log C/<-0.82 

16.63 

16.54<log N c i V ^16.80 

19.6T0.1 

0 qq+0.59 
^•^-0.47 

21 

9 

-0.84 

-0.96<log U<-0.78 

16.60 

16.52<log N c i V ^16.77 

18.6T0.1 

0 228~’~°‘° 59 
u.zzo _ 0 Q47 

+22 

5 

0.08 

-0.02<log t/<0.12 

16.77 

16.70sSlog N c i V ^17.03 

20.1T0.1 

72 o+ 18 - 7 
/z - u -14.8 

22 

7 

0.06 

-0.04<log (7^0.10 

16.74 

16.67^1og N c i V ^17.04 

19.2T0.1 

9 05”*" 2 ' 35 

1.85 

22 

9 

0.08 

-0.02<log C/^0.14 

16.72 

16.62$log N c iv<16.95 

18.1T0.1 

0 720^”°' 187 
u - / zu _ 0. 148 

+23 

5 

1.12 

l.OO^log CAg 1.16 

16.66 

16.59^1og N c iv<16.90 

19.6T0.1 

228t®?* 

23 

7 

1.08 

0.98^1og t/^1.14 

16.68 

16.57sSlog N c iv<17.00 

18.6T0.1 

22.8^® 

23 

9 

1.12 

1.00<log U ^1.16 

16.61 

16.54^1og N c iv<16.81 

17.6T0.1 

r\ qq + 0.58 
^•^0-0.47 


* These outflow rates are much greater than the accretion rate and are therefore ruled out. 
' These Nh,Rh grid-points are ruled out due to the large outflow rates indicated by *. 


e.g. in lReimers et al] (120051 '). where the flux at a rest-frame 
wavelength of 335 A changed by a factor of ~3 in a rest- 
frame timescale of ~0.65 years. Since the EUV is the main 
contributor to outflow ionization, this indicates our results 
are plausible. 

As indicated in Column 7 of Table [7] the log R values 
do not change significantly between epochs for given combi¬ 
nations of log Nh and log ?ir. The dependence of logR on 
log Nh is indicated in Fig. 1101 It is appropriate to consider 
how the predicted absorber distan ces compare to the size 
of the BLR. Using the method of iKaspi et ali (l2005lb the 
BLR radius Rblr can be estimated using Rblr= AL s x 10 
lt-days, where L is AL>(1450 A)/10 44 (erg s _1 ) and A and B 
are constants, with a set of values found using t he F I TEXY 
method and the mean Balmer line’s time lag in lKaspi et al.l 
J2005ll to be A=2.12 and B=0.496. This set of values is used 
as the index B is consistent with subsequent studies, e.g. 
[Bentz et alj (|2006h which found it to be ~0.5 This is consis- 
tent with the fact that UocQh/R 2 , leading to RocL 0,5 . The 
procedure provides a value of .Rbl.r~7.6x 10 17 cm (a few 


hundred light-days). Comparing this with the values found 
for R in Table [7] suggests that the outflowing gas is out¬ 
side the BLR in 8 out of 9 cases, the exception being the 
grid-point log (Nh / cm _2 )=23, log (uh / cm _3 )=9. 

The mass outflow rates calculated range from tenths 
of a solar mass to hundreds of solar masses per year. The 
estimated mass outflow rates, which depend upon the prod¬ 
uct RNh, are in the range indicated in Table [3 Given the 
uncertainties in the mass outflow rate, they do not change 
significantly between epochs. It is possible to calculate the 
accretion rate at Eddington luminosity in this object and 
hence compare it to the actual accretion rate and the mass 
outflow rate. This requires an estimation of the black hole 
mass, which can be calculated using the following 

M = (8 ) 

Gr 

where v is the orbital velocity of the BLR gas, G is the grav¬ 
itational constant and / is the virial coeffi cient which has 
been estimated to be ~5.25 (I Woo et all201Ch . The value of v 
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Figure 10. Dependence of log R on log Nh for constant mass 
outflow rates. Values are from Table [3 

is estimated as v~4100 km s -1 using the width of the C IV 
broad emission line. This gives a value of Mbh ~5.1x 10 9 
Mg. This was checked for consistency by compa ring with 
the value obtained using Equation 3 in I Vestergaardl (1200411 . 
which provides a very similar value of ~5.3xl0 9 Mg. The 
Eddington luminosity can then be found using the relation 
LEdd=W i61 Mg erg s -1 , where Ms is black hole mass in 
units of 10 8 Mq, giving Lsdd= 6.4xl0 4 ' erg s -1 . The Ed¬ 
dington accretion rate MEdd=LEdd/i ?c 2 implies Mndd^H 3 
Mq yr -1 . 

The actual accretion rate can be calculated from 
the bolometric luminosity, which can be found using 
the lumi nosity at 145 0A by applying Li, 0 ;~4.36AL> (at 
1450A) (1 Warner et alj 12004 1. For a bolometric luminosity 
of Tt, 0 i~9.5x 10 46 erg s~, the resultant value of the ac¬ 
cretion rate in terms of that at Eddington is ~0.15 Mndd 
(~17.0 Mgyr -1 ), meaning that this object’s accretion rate 
is a substantial portion of the value at the Eddington limit. 
This allows the ruling out of two combinations of hydrogen 
column density and hydrogen number density (log (Nh / 
cm -2 )=22, log ( n,H / cm -3 )=5 and log (Nh / cm -2 )=23, 
log (n,H / cm -3 )=5) due to their excessively high outflow 
rates compared to the accretion rate, as indicated in Table[7| 
Comparing the remaining allowed values suggests the mass 
outflow rate is approximately between 1 per cent and 200 per 
cent of the mass accretion rate, corresponding to a kinetic lu¬ 
minosity range of~3.3xl0 42 to6.6xl0 44 ergs -1 . The upper 
limit of this range gives a kinetic energy to luminosity ratio 
of ~0.7 per cent, high enough to have a significant impact 


on th e evolution of the surrounding galaxy dHopkins fo Elvii 
l20ldl . however definite conclusions regarding the outflow’s 
importance in this regard are impossible given the weakness 
of the constraint. 


5 DISCUSSION 

5.1 Are changes in the outflow geometry 
responsible for variability? 

Changes in covering fraction within the context of the PPC 
model are one possibility for the dramatic variability seen in 
this quasar’s absorption lines. From Section 3.4 (Fig. [4]) it 
is clear that the deepest parts of the blue absorber 1 trough 
in epoch 2 are accurately described by a saturated absorber 
whose depth is determined by the fraction covering the emis¬ 
sion source. However, Section 3.5 (Fig. (5j) indicates that an 
inhomogeneous absorber modelled over the same velocity 
range is also an excellent predictor of the same trough. Fur¬ 
ther insight can be gained from the Gaussian model profiles, 
specifically absorbers 1 and 2 which are well-defined and 
hence can be studied using both PPC and inhomogeneous 
models. From Tabic [o] it appears that these features are not 
well described by an inhomogeneous absorber, since, given 
the maximum power-law index of 10, many of the calculated 
values are unconstrained. Large power indices produce line- 
of-sight geometries which resemble a step-function, which 
is equivalent to a PPC absorber. The difference between 
the power-law index which matches the deepest part of Ab¬ 
sorber 1 in epoch 2 (<z=3.3+q' 4 ) and the value in Table 0 
is also much greater than the uncertainty. This leads to the 
conclusion that an inhomogeneous absorber is probably not 
the predominant influence on absorption line profile shape. 

If it is assumed that absorption takes place within the 
context of PPC, it is possible to determine how the cover¬ 
ing fraction has changed between epochs using the Gaus¬ 
sian profile models. From Table [4] is is clear that there is 
no significant change in C for absorber 1 between epochs 
(changes being within la for <7^0.042). There is also no dis¬ 
cernible change in the covering fraction of absorber 2. Con¬ 
clusions regarding covering fraction variations in absorber 
2 are weak for two reasons. First, the epoch 2 doublet is 
entirely unsaturated and therefore the covering fraction es¬ 
timation is completely unconstrained. Secondly, when com¬ 
paring the epoch 1 value (C=0.187±0.017) with the epoch 3 
value (C=0.275±0.165), the large uncertainty in the latter 
allows a significant range which prevents any strong conclu¬ 
sion regarding the presence of covering fraction variations. 
Given the lack of significant variability in either the strongly 
constrained absorber 1 or the weakly constrained absorber 
2 , we conclude that there is no evidence for covering fraction 
being the dominant mechanism behind the variability. 

5.2 Evidence for ionization changes driving 
variability 

Ionization fraction changes within the outflowing gas can be 
generated either by intrinsic changes in the ionizing radia¬ 
tion emitted from the central ionizing continuum source or 
shielding of this radiation by gas in the region between the 
accretion disk and the BAL gas. There are several indicators 
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that suggest ionization fraction changes are implicated in the 
absorption line variability in this object, the strongest being 
coordinated changes seen ac ross la r ge v elo city se parations. 
A similar result was noted in lTrevese et all (12013h regarding 
the quasar APM 08279+5255, where they observed coordi¬ 
nated variability in two Gaussian-modelled components of 
a C IV BAL separated by ~5000 km s -1 . Large velocity 
separations between absorption troughs imply that the ab¬ 
sorption arises in physically distinct regions, hence making 
covering fraction changes an unlikely explanation due to the 
level of coordination requ ired across such large separations 
(e.g. iFiliz Ak et al.l ll2013T )~). 

This phenomenon is especially evident in the C IV ab¬ 
sorption troughs (see Fig. H. The visibly saturated absorp¬ 
tion trough centred at 4000 km s -1 undergoes weakening 
between epochs 1 and 2 before strengthening again in epoch 
3, which is a pattern also seen in the features at ~—9500 km 
s' 1 and ~—18 500 km s _1 . Since there is a lot of blending 
of the Si IV components, the same effect is not as obvious in 
this ion. However, it is instructive to look at the behaviour 
of the modelled Gaussian Si IV doublets in Fig. [3] where the 
same effect is apparent. 

No part of the power-law continuum contributing to this 
spectrum is observable short-ward of Ly-a, however it is rea¬ 
sonable to exp ect that it extend s in to the EUV portion of 
the spectrum dKriss et alJll999h . which is largely responsi¬ 
ble for the flux which controls the ionization fractions of the 
plasma in the accretion disk’s environment. Therefore vari¬ 
ability in the observed continuum could indicate variability 
in the ionizing radiation from the accretion disk. There is 
evidence for an excess in the blue part of the spectrum dur¬ 
ing epoch 2 (first WHT observation) when compared to the 
prior and subsequent observations, as indicated in Figs 1111 
and ED which could indicate an increase in the ionizing con¬ 
tinuum flux at this epoch. At the 3rd epoch the blue excess 
appears to disappear, with the flux in this range appearing 
to follow the continuum from epoch 1 once again. 

Evidence from the CLOUDY simulations also supports an 
ionizing continuum change driving the variability. Given the 
decrease in the Si IV column density between epochs 1 and 
2 , an increase in ionization parameter must occur between 
these observations (See Figures [3] and [3]), with an estimated 
log U change of ~0.5 to ~0.6 dex. An increase in log U 
agrees with the observed steepening of the continuum at the 
blue end of the spectrum between epochs 1 and 2, assuming 
the continuum extends to ionizing energies. There is a large 
overlap in predicted Si IV column densities between epochs 
2 and 3, however the epoch 3 range lies at the higher end 
of the epoch 2 range, with the best epoch 3 value being 
0.3 dex above the best epoch 2 value. A log U change of 
—0.2 dex is needed to drive this increase in Si IV column 
density between epochs 2 and 3. A decrease in ionization 
parameter between these two epochs is supported by the 
apparent disappearance of the aforementioned blue excess 
in the continuum. 

5.3 Fourth Observation 

The quasar SDSS J1138+3517 was observed again by the 
WHT in May 2014 and is included in Fig. [3] In contrast 
to the earlier WHT observations, the Si IV components 
are largely unchanged from the previous epoch’s spectrum. 



4000 5000 6000 7000 


Observed Wavelength (A) 

Figure 11. Upper panel: Epoch 1 observation (black line) and 
epoch 2 observation (red line) with power-law continuum (dashed 
black line) fitted to the epoch 1 spectrum. Vertical dotted lines 
indicate the boundaries of where the blue excess is apparent. The 
thick black horizontal line indicates the dichroic overlap region 
in the WHT observation. Lower panel: Same as upper panel but 
with red line indicating the epoch 3 observation. 


The deep C iv features centred at ~—4500 and ~—9500 
km s -1 do not show changes comparable to previous vari¬ 
ability when going from epoch 3 to epoch 4. However the 
weaker feature at ~—18 500 km s _1 shows a significant 
strengthening, with its profile becoming comparable to that 
of the epoch 1 observation. The absorption line centred at 
~—13000 km s , which is very conspicuous at epoch 1 but 
almost unidentifiable at epochs 2 and 3, shows a slight in¬ 
crease in strength, however the significance of this change is 
only 1.08<r. 

Although these changes from epoch 3 to epoch 4 are in 
the weakest C IV components and therefore do not signifi¬ 
cantly change the minimum column density estimate made 
in epoch 3, the localisation in velocity space of this strength¬ 
ening to a range of ~6000 km s -1 is interesting. If the in¬ 
crease in absorption is only considered to be significant in 
the stronger of the two features, then the region over which 
the change occurs is reduced to ~4600 km s -1 . A single 
varying trough allows the possibility of a change in the cov¬ 
ering fraction of the absorber to explain the phenomenon. 
This trough is also interesting in that it does not show as 
great a change between epoch 2 and 3 as the feature centred 
at —9500 km s -1 . This could be due to either (a) the line 
being the result of a moderately saturated outflow whose 
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Rest —frame Wavelength (A) 

1250 1300 1350 1400 



Observed Wavelength (A) 


Figure 12. Epoch 1 observation (black), epoch 2 observation 
(red) and epoch 3 observation (green), normalised to the epoch 
1 continuum shown in Figure ITT! over the blue excess range, with 
dotted lines representing the errors. The epoch 2 spectrum ap¬ 
pears to lie almost everywhere above the epoch 1 continuum. 

peak trough depth describes the peak covering fraction, or 
(b) the outflow is located in a region where it is relatively 
shielded from the ionizing flux variations (the large velocity 
separation from the other troughs suggests it is in a physi¬ 
cally distinct location). 


6 SUMMARY 

This following provides a summary of the conclusions that 
can be drawn from this investigation. 

(i) When c omparing to large-sample s tatistical studies o f 
BALQSOs dFiliz Ak et al.l 1 20131 . 1 20141 : IWildv et al.l I2014I) . 
BALs in SDSS J1138+3517 are exceptionally variable, espe¬ 
cially in the 360 day rest-frame time interval between epoch 
1 and epoc h 2. When measuring BAL variability using the 
method of iFiliz Ak et all d2013l . 120141 ). the most variable 
C iv BAL in this object undergoe s a change of AEW= — 13.8 
A. When compared to Fig. 29 of IFiliz Ak et al.l (120131 ). this 
change in A EW is far outside the range predicted by their 
random-walk model at an equivalent time interval. This 
could indicate that SDSS J1138+3517 belongs to a rare class 
of highly variable BALQSOs whose variability mechanism is 
not the s ame as t hat driving the random-walk behaviour ob¬ 
served in iFiliz Ak et al.1 (|2013l ). 


(ii) . When applying Gaussian fitting to model the Si iv out¬ 
flows in Section 3.3, it is found that the components of in¬ 
dividual doublets rarely exhibit the expected 2:1 blue to 
red optical depth ratio and also do not reach zero intensity 
at line centre. This suggests saturation of the blue compo¬ 
nent combined with spatially dependent line-of-sight optical 
depth (either inhomogeneous coverage or pure partial cov¬ 
erage). 

(iii) . Although both inhomogeneous power-law and PPG 
models successfully predict the central velocity range of the 
deepest saturated blue trough in epoch 2, modelling the ab¬ 
sorption as Gaussians over the entirety of the velocity range 
suggests a PPC model is more successful. This is due to the 
fact that the best values of the power-law index suggest a 
very steep slope, therefore approximating a step function 
which is the mathematical equivalent of PPC. The fact that 
both models can be successfully applied over a limited range 
in absorption line velocity is something that should be noted 
in future investigations which attempt to distinguish these 
models. 

(iv) . There is no evidence for a change of covering fraction 
across the epochs when calculating C from the (Gaussian 
modelled) doublets representing the two deepest absorbers, 
suggesting line-of-sight covering changes are not the domi¬ 
nant mechanism driving the absorption line variability. How¬ 
ever the weak C IV feature centred at 18 500 km s^ 1 may 
have undergone a covering fraction change between epochs 
3 and 4. 

(v) . There is evidence that ionization changes play a signif¬ 
icant role in the variability, most notably from the coordi¬ 
nated changes in both Si IV and C iv absorption troughs 
separated widely in velocity. There is also a hint that the 
ionizing flux may have increased in epoch 2, given the ap¬ 
parent steepening of the continuum in the blue part of 
the spectrum. This conclusion would be supported if AGN 
have steeper a ox a t higher luminosities as indicated in 
IVignali et alj (120031) . The photoionization simulations are 
consistent with this scenario, since the column densities of 
both ions decrease with increasing ionization parameter in 
the log U range spanned by the outflow. The best values for 
the changes in ionization parameter are within the scope of 
quasar EUV variability. 

(vi) . The location of the inner face of the outflowing gas is 
not well constrained, being located between approximately 
0.1 and 200 parsecs from the continuum source. This range 
places it somewhere in a region spanning the BLR radius 
to narrow line region distances. However, the lowest veloc¬ 
ity (and deepest) C IV trough must be absorbing both BLR 
and continuum emission, meaning that this outflow is lo¬ 
cated outside the BLR. Given the width of the absorption 
components, which are substantially less than BLR widths, 
and the value of Rblr, it seems likely that the inner face is 
at least ~1 parsec from the continuum source. 

(vii) . The magnitude of the mass outflow rate could be a sub¬ 
stantial fraction of or even exceed the mass accretion rate. 
However due to the uncertainties on R, strong constraints on 
the true value of the mass outflow rate are not achieved. The 
derived kinetic luminosity range of ~3.3xl0 42 to 6.6xl0 44 
erg s _1 indicates that feedback effects such as quenching of 
star formation may be possible if the true value is near the 
upper limit of this range. 
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